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ABSTRACT 

Annihilation of dark matter can result in the production of stable Standard Model particles including 
electrons and positrons that, in the presence of magnetic fields, lose energy via synchrotron radiation, 
observable as radio emission. Galaxy clusters are excellent targets to search for or to constrain the 
rate of dark matter annihilation, as they are both massive and dark matter dominated. In this study, 
we place limits on dark matter annihilation in a sample of nearby clusters using upper limits on the 
diffuse radio emission, low levels of observed diffuse emission, or detections of radio mini-haloes. We 
find that the strongest limits on the annihilation cross section are better than limits derived from the 
non-detection of clusters in the gamma-ray band by a factor of ~ 3 or more when the same annihilation 
channel and subtructure model, but different best-case clusters, are compared. The limits on the cross 
section depend on the assumed amount of substructure, varying by as much as 2 orders of magnitude 
for increasingly optimistic substructure models as compared t o a smooth N FW profile. In our most 
optimistic case, using the results of the Phoenix Project (jGao et al.ll2Q12bl ). we find that the derived 
limits reach below the thermal relic cross section of 3 x 10 -26 cm 3 s _1 for dark matter masses as large 
as 400 GeV, for the bb annihilation channel. We discuss uncertainties due to the limited available data 
on the magnetic field structure of individual clusters. We also report the discovery of diffuse radio 
emission from the central 30-40 kpc regions of the groups M49 and NGC4636. 

Subject headings: dark matter - galaxies: clusters: individual: M49, NGC4636 - galaxies: clusters: 
intracluster medium - radiation mechanisms: nonthermal - radio continuum: gen- 
eral 

1. INTRODUCTION 

Clusters of galaxies are the most massive virialized ob- 
jects in the universe. About 80% of the mass of clusters 
is comprised of dark matter, making them good candi- 
dates for astrophysical searches for a signature from par- 
ticle dark matter. Among the best motivated particle 
candidates for dark matter are weakly interacting mas- 
sive particles, or "WIMPs". WIMPs can self-annihilate 
to Standard Model particles, including electrons and 
positrons. Example particle theories that predict WIMP 
dark matter are supersymmetric e xtensions to the Stan- 
dard Model (for a review, see e.g. JJungman et aTI [l996) 
and universal extra dime nsions (for a review, see e.g., 
IHooper Profumol [2007h . The products of WIMP an- 
nihilations generically yield a broad spectrum of electro- 
magnetic emission in cluster environments that is poten- 
tially observable acros s a wide range of frequenc ies, from 
radio to gamma rays (Colafr ancesco et aI1l2006f ). 

Recent searches for dark matter annihilation have 
largely focused on indirect detection through gamma-ray 
emission. Gamma rays are primarily produced in WIMP 
annihilation by two classes of processes: 

(i) in the prompt emission from the two-photon de- 
cay of neutral pions produced in the hadronization of 
strongly-interacting particles produced in the WIMP an- 
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nihilation, or from radiative processes such as internal 
bremsstrahlung, and 

(ii) in the secondary emission when high-energy elec- 
trons and positrons (either directly produced or result- 
ing from the decays of particles produced in the anni- 
hilation process, e.g. charged pions) up-scatter ambient 
photons to gamma-ray energies, a process known as In- 
verse Compton (IC) scattering. 

An additional secondary gamma-ray emission pro- 
cess, typically subdominant, is associated with secondary 
bremsstrahlung off of ambient gas of the same elec- 
trons and positrons resp onsible for the IC emission 
(jColafrancesco et al.ll2006l ). 

Galaxy clusters have not yet been detected in gamma 
rays (see recen tly, e.g., iMacfas- Ramirez et al I [20121 : 
lHan et aTl l2012f ). Upper limits on the gamma ray 
emission have been used to place constraints on the 
dark matter annihilation cross section. Previous stud- 
ies that placed constraints on dark matter annihila- 
tion or decay in clusters have focused on searching 
for gamma-ray emission using the Large Area Tele- 
scope (LAT) onboard the F ermi Gamma Ray Space 
Telescope at the GeV scale (Ackermann et al. 2 010al: 
i Dugger et all l2010[ lAndo fc N agai 2012; iHuang et all 
120121 : iNezri et al.l 120121 : lHan et all 120121 ). and using 
ground-based Cher e nkov telescopes, inclu ding H.E.S.S. 
(lAbramowski et al.l l2Q12h and MAGIC (lAleksic et al.l 
2010) at the TeV scale. 

The leptons that IC-upscatter photons to GeV and 
TeV energies also lose energy via synchrotron radiation in 
the presence of magnetic fields, producing emission that 
is potentially observable at radio frequencies. Diffuse 
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radio emission from the intracluster medium (ICM) in 
the form of approximately spherically symmetric haloes 
and mini-haloes has been observed from more than 40 
clusters to date, indicating the existence of both a rela- 
tivists population of electrons/positrons and large scale 
mag netic fields at the ~ fiG level (for a review see 
e.g., Fer etti et al.l 120121 ). Giant radio haloes are typi- 
cally ~ 1 Mpc in size, tend to follow the distribution 
of the thermal X-ray emission and are therefore roughly 
spherical, and are characterized by steep energy spectr a 
and low surface brightness (e.g.. iThierbach et al.l 12003). 
The first discovered, and perhaps best studied mem- 
ber of this class is th e halo of the Coma Cluster (see 
e.g. [L arge et al. 1959; Willson 1970; iKim et al.l [199(1 



iKronberg et al.l 120071 : iBrown fc Rudnickl 1201 ID . Mini- 
haloes are similar in shape to haloes but are on the order 
of hundreds of kiloparsecs in size and are found in clusters 
with cool cores, which are characterized primarily by a 
very short central cooling time (<C 1 Gyr), and addition- 
ally show a sharp increase in X-ray surface brightness and 
a corresponding d ecrease in temperat ure at the center of 
the cluster (e.g., Hud son et al.ll20To[ ). The prototype of 
this class is the m ini-halo found in the Perseus Cluster 
([Burns et al.lll992f ). 

The origin of these relativistic particles responsible for 
the diffuse radio emission is, however, at present unclear. 
They perhaps were injected into the ICM by AGN or 
supernova activity within cluster galaxies and then reac- 
celer ated by turbulence in the ICM caused by mergers 
(e.g.. iBrunetti et al.l 120011: iPetrosianl 120011 : iFuiita et all 
2003; iCassano fc Brunettill2005[ ). Thev could also be the 
products of cosm ic ray hadron collisions with particles 
in th e ICM (e.g.. iDen nison 1980; Blasi & Colafran cescol 
1999; iPfrommer fc EnBlin 2001)7 This scenario requires 
that magnetic fields be strong, at least comparable 
to the strengths inferred from Faraday Rotation Mea- 
sures (RMs), so that the gamma ra y fluxes from cosmic 
rays do not exceed current lim its (jJeltema fc Profumol 
120111 : lAckermann et al.l l2Q10bh . Finally, the relativis- 
tic electrons and positrons could also be the result 
of dark matter annihila tion (Colafra ncesco et al.l 2006: 
iPerez-Torres et al.ll2009l ). 

Most clusters show no or only lo w levels of de t ectable 
diffuse radio emission. However, IBrown eTaD (j2011bh 
stacked a sample of 105 clusters with no detected radio 
haloes and found a 6a detection of radio emission, which 
suggests that even clusters with no currently observed 
radio emission host, at some level, both relativistic pop- 
ulations of electrons and positrons and large-scale, jaG 
magnetic fields. The relativistic electrons and positrons 
potentially produced by dark matter annihilation in clus- 
ters would synchrotron radiate away their energy in the 
presence of cluster-strength magnetic fields. We can thus 
use the radio upper limits or low levels of radio emission 
from clusters to constrain the rate of dark matter annihi- 
lation. In this paper, we place limits on the dark matter 
annihilation cross section using primarily upper limits on 
or low levels of the diffuse radio emission from a selec- 
tion of nearby clusters and also using some clusters with 
detected mini-haloes. Our study greatly improves con- 
traints on dark matter with radio observations over pre- 
vious studies by carefully selecting optimal targets with 
little or no observed diffuse radio emission, as opposed 
to clusters with bright radio halos that have been con- 



sidered so far (see e.g. Colafr ancesco et aT1l2006[ ). 

The paper is organized as follows. In Sec. [2j we de- 
scribe our models for the dark matter and magnetic field 
profiles in clusters. In Sec. [3] we present the sample of 
clusters and radio data we employ for our analysis, re- 
port on the observation of diffuse emission from M49 
and NGC4636, and discuss the impact of studies of intr- 
acluster magnetic fields on our sample. Sec. |4] presents 
limits on the dark matter annihilation cross section and 
discusses the implications and caveats of our results, in- 
cluding comparisons to previous studies and effects of un- 
certainties in the assumed dark matter profiles and mag- 
netic field structure. We conclude in Sec. Through- 
out this paper, we assume a ACDM cosmology with 
H = 100k km s -1 Mpc -1 , h 0.70, Q m = 0.27, and 
Q A = 0.73. 

2. DARK MATTER AND MAGNETIC FIELD MODELING 

The spectral flux density due to dark matter annihila- 
tion is: 
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where (av) is the thermally- averaged zero-temperature 
dark matter annihilation cross section times velocity, 
m x is the dark matter particle mass, and EdN v /dE is 
the synchrotron energy spectrum per frequency per dark 
matter annihilation event. The observable flux density 
S v is typically measured in Jansky (Jy) in the radio. J 
is the line-of-sight integral of the dark matter density 
squared, integrated over the angular size of the emission 
region AQ: 
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We take the size of the emission region to match that of 
the available radio data, either the size of the observed 
diffuse emission or the size used to place upper limits. 
We use a number of different models for the dark matter 
density profile, described in the following sections. 

2.1. NFW Density Profile 

In our analysis we use four models for the dark mat- 
ter density profile. The most conservative, with no 
substructure, is the Navarr o- Frenk- White (NFW) profile 
(jNavarro et aLlll996L fl997h : 



PNFw(r) 
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where the central density p s and the scale radius r s are 
determined by ob s ervatio ns. Following the derivation in 
lAckermann et all ([201 0a), we use the scaling relation- 
ship derived from X-ray observations of clusters from 
iBuote et al.l ([20071 ) to determine r s and p s from the virial 
mass and radius, M vir and r V i r : 
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where the concentration 

Gvir — ^vir/^s' The virial over- 
density A vir « 98 is defin ed with respect to the critical 
density of the universe (|Brvan fc Normanl Il998h , with 
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M vir = ^- A vir p c r^ ir . We determine M vir and r V i r from 
M500, where the overdensity is 500 times the critical 
density, obtained from X-ray observations ([Chen et al.l 
|2007[ ). We correct M500 for our current cosmology and 
reduce it by the gas fraction f gas , also reported in 
I Chen et al.l ([20071 ). so that we are left with only the mass 
of the dark matter in the cluster: M500 — > M5oo(l — 
fgas)h. We use the corrected M500 to determine r$oo 
and use these parameters to determine M^ r , r V i ri and 
r, s using the follow i ng eq uations, derived in Appendix C 
of lHu fc Kravtio^l (|2003f ): 
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A bm f(r s /r vir ) = A vir f(r 3 /r 50 o), 
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f(x) = x 3 [ln(l + X' 1 ) - (1 + a?)" 1 ]. (7) 

We u se M V i r , r V i r , and r s to solve for p s ()Hu fc Kravtsovl 
G2QQ3): 



^f(r s /r vir ) 

2.2. Effects of Substructure 

In cold dark matter (CDM) cosmologies, dark mat- 
ter haloes are structured hierarchically. Clusters are 
observed to host subhaloes down to at least the scale 
of dwarf galaxies, approximately 10 T M©, and are pre- 
dicte d to contain sub structure down to ~ 10~ 6 M @ 
(see lGreen etHI (j2005f ). but also lProfumo et all ([2006( 1). 
Since the J factor for dark matter annihilation is propor- 
tional to the density squared, the amount of substructure 
is critical, and can increase J by one or more orders of 
magnitude compared to a smooth NFW density profile. 
In clusters, tidal stripping tends to destroy subhaloes 
near the cluster center, so substruct ure is preferentiall y 
found towards the cluster outskirts ([Gao et al.l 120 12b). 
Generally the density profile of subhaloes is more ex- 
tende d than the smooth central halo, as seen in simula- 
tions (jNagai fc Kravtsovll2QQ5l iGao et al.ll2Q12bh . 

We consider here multiple s ubstructure models . Two o f 
these models are also used in lAckermann et al.l (|20 10a): 
a conservative model, with a subhalo cutoff mass of 
10 7 Mq and a fraction of the total halo mass in sub- 
structure, f Sl of 10%, and an optimistic model, with a 
cutoff mass of 10~ 6 M© and ,f, s = 0.2. We use the formal- 
ism in Colafr ancesco et aTl (|2006l ) to determine the dark 
matter density profiles pcon and popt- In the notation 
of this framework, the conservative setup corresponds to 
a substructure contrast factor of A 2 = 1.3 x 10 5 , and the 
optimistic setup corresponds to A 2 = 3.0 x 10 5 . 

We also consider a third model based on the results 
of the Phoenix Project, which is a series of dark mat- 
ter simul ations following the evolution of cluster-sized 
haloes ([Gao et al . 20 12a. 3). These simulations adopt 
the cosmological para meters of the Millenium Simulation 
([Springel et al.ll2005f ). which are now known to be incon- 
sistent with current cosmological parameters. However, 
the main difference is the value of ag, which affects the 
number of clusters in a universe-sized simulation but not 
the properties of individual cluster halos. This model 



adopts a cutoff subhalo mass of 10 -6 M , resulting in a 
subhalo mass fraction of approximately 27%, and is more 
extended than either of the previous models. 

The Phoenix simulations use M200 and r2oo to express 
their results in ter ms of a boost fact or and substructure 
surface brightness (Gao et al. 2 012at ): 
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where L ma i n is the luminosity of the smooth halo, which 
is well-described by an NFW profile. We use Eqns. (|5]) 
and (|6]) to solve for M200 and r2oo> defined through the 
overdensity A200 = 200p c , from M500 and rsoo- We then 
translate Eqn. (fTOl) to a J- factor "surface brightness" , 
following Ha n et aD {2012). Thus the resulting J factor 
for the Phoenix simulations is: 



Jphx = Jnfw + Jsub 



where 
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Electron/Positron Signal from Dark Matter 
Annihilation 



In order to characterize the synchrotron spectrum, 
we need the equilibrium electron and positron spectra, 
which result by solving the full diffusion equation: 
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where dn e /dE is the equilibrium electron or positron 
density spectrum, Q(E,x) is the source term, D(E,x) 
is the spatial diffusion coefficient, and bi oss (E,x) is the 
energy loss term described below. In clusters of galaxies, 
the timescale on which electrons radiate is much shorter 
than the spatial diffusion timescale, so we can neglect 
the time dependence on the l.h.s. as well as the spatial 
depe ndence on the r.h.s. of Eq . (fT3|) (see e.g., Appendix 
A offColafran cesco et al. 2006 for a discussion of the role 
of diffusion in clusters). The source term is proportional 
to the injected spectrum of electrons/positrons per dark 
matter annihilation. The expression for the equilibrium 
density spectrum is: 



dn e 
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We use the DMFIT package ([Jeltema fc Profumol 
2008) (in turn deri ved from the DarkSUSY package 
([Gondolo et al.ll2004l )) to calculate the electron/positron 
injection spectra per dark matter annihilation. The en- 
ergy loss term, in the low redshift limit, is the sum of 
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synchrotron, IC, bremsstrahlung, and Coulomb losses: 



bioss(E) 



blC + hrem + b c 



0.0254 
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+ 1.51n(0.36 + log(7/n)) 
+ 6.13(1 + log(7/ra)/0.75) 

The energy loss term bi oss (E) has units of 1 x 10 -16 GeV 
s _1 , where n is the average thermal electron density, 
« 1 x 10 -3 cm -3 for all clusters. For GeV electrons 
and positrons, synchrotron and IC losses dominate; when 
B > Be mb ~ 3 /iG, synchrotron losses dominate over 
IC losses. 

2.4. Synchrotron Emission 

In clusters, relativistic electrons radiate their energy 
via synchrotron emission in the presence of a mag- 
netic field B(r), and live in a background plasma 
with electron density n(r) and plasma frequency v v = 
8890[n(r)/l cm -3 ] 1 / 2 Hz. The power per frequency of 
emitted synchrotron radiation for a single electron (or 
positron) with energ y E = 7m e c 2 , averaged over all in- 
coming directions is (|Longairll2011[ ): 
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where ro = e 2 /(mc 2 ) is the classical electron radius, 
is the pitch angle, v$ = eBj(2i\mc) is the nonrelativistic 
gyrofrequency. The quantities x and F are defined as 
follows: 



2v 



x = 



3^o7 2 

« 1.25,s 1/3 exp(- ; s)[648 



3/2 



(17) 



(18) 



2il/12 



■A 



where if 5 / 3 (£) is the modified Bessel function of order 
5/3. 

The synchrotron energy spectrum per frequency, 
EdN u /dE, given populations of electrons and positrons 
each with an equilibrium density spectrum dn e /dE e is: 
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This energy spectrum is inserted into Eq. (pQ) to find the 
limits on the annihilation cross section. 

2.5. Magnetic Field Model 

Synchrotron emission depends strongly on the mag- 
netic field in the region of interest. From MHD simu- 
lations and from clusters with multiple Faraday RMs, 
the magnetic field appears to follo w the gas density 
(iMurgia et al.l 120041 : IBonafede etHI 120101 : iVacca et all 
120121 )! The gas density is typically fit with a f3 model 



(|Cavaliere fc Fusco-Femiano| [r976): we adjust the radial 
dependence with a free parameter n: 



B(r) = B 
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where Bq is the central magnetic field value and r c is the 
core radius. The parameters j3 and r c are fit using X- 
ray data, and Bo and n are typically modeled based on 
a combination of sim ulations and RMs (e.g., for Coma, 
IBonafede et~aT1l2010[ ). 

We weight the magnetic field distribution by the dark 
matter density profile to yield a single parameter, an 
effective magnetic field B e ff that is used in our syn- 
chrotron spectrum calculations: 
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where is the radius of the considered emission region, 
either that of the observed diffuse emission region or the 
region used to place an upper limit. This weighting of 
the magnetic field yields a better estimate of the field 
in the regions where most of the synchrotron emission 
is originating from, i.e., regions with high dark matter 
densities. 

3. CLUSTER SAMPLE: RADIO AND MAGNETIC FIELD 

DATA 

3.1. Selection of Clusters from Radio Data 

We choose our sample based on radio, X-ray, and mag- 
netic field data available from the literature. The clus- 
ters that produce the best limits on dark matter an- 
nihilation are nearby and do not host observable radio 
haloes, or host only low levels of central diffuse emis- 
sion or mini-haloes. Unfortunately, there are very few 
published upper limits on the diffuse ra dio emission in 
clusters. iRudnick fc Lemmermanl (|2009l ) published the 
most comprehensive list of upper limits (96% confidence 
limit) on Mpc-sized radio emission for bright X-ray clus- 
ters with redshifts between 0.03 and 0.3 at 327 MHz with 
the Westerbork Northern Sky Survey. Upper limits on 
Mpc-sized e mission from clusters with 0.2 < z < 0.4 were 
published in lVenturi et al.l (|2008[ ) using the Giant Metre- 
wave Radio Telescope at 610 MHz. Since the astrophys- 
ical J factor decreases with increasing distance, nearby 
clusters are the best candidates for dark matter detec- 
tion; we therefore restri ct our sample to z < 0.1 . Only 
a handful of clusters in IRudnick fc Lemmermanl ([2009) 
have redshifts of less than 0.1. Additionally, to facil- 
itate comparisons, our sample of clusters largely over- 
laps with those considered for gamma-ray a nalysis (e.g., 
lAckermann et al.ll2010al iHuang et al.ll2012f ). 

We take f3 model parameters, cluster masses, and gas 
fractions from the HIFLUGCS catalo g of nearby, mas - 
sive, bright X-ray galaxy clusters ([Chen et all 12007). 
Four clusters with publis h ed rad io halo upper limits from 
IRudnick fc Lemmermanl (|2009f ) are also part of the HI- 
FLUGCS catalog. 

We also include three additional clusters in HI- 
FLUGCS that have detected mini-haloes: Perseus, Ophi- 
uchus, and A2029. Despite the detected radio flux from 
these clusters, they are strong candidates because they 
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Figure 1. Greyscale image of M49 from NVSS, 45" resolution, 
on a scale of 0-100 mJy/45" beam. Contours show the residual 
diffuse emission at a resolution of 135", as shown in the lower left, 
after subtraction of the compact emission. Contour levels are at 2 
mJy/135" beam x (1,2,3,...). 

host cool cores, and thus have higher inferred central 
magnetic field strengths, and are also nearby and heav- 
ily dark matter dominated. Perseus and Ophiuchus es- 
pecially are considered to be among the best candidates 
for dark matter searches ( Jel tema et al.l 120091 ). Coma is 
also included in ou r sample, as it h as both a very well- 
studied radio halo (iKim et al.|[l990l ) and magnetic field 
(jBonafede et al.ll2010f ). and is generally used as a rep- 
resentative massive, merging cluster. It is also the only 
cluster that has been used previously to place constraints 
on the dark matter conten t in clusters using radio data 
(jColafrancesco et^ al. 2006) and is frequently considered 
in analyses of dark matter annihilation in clust ers, espe- 
cially in gamma rays (e.g., most recently in lHan et al.l 
|2012[ ). We show that it produces comparably poor lim- 
its on the dark matter annihilation cross section due to 
the large observed radio emission, for a smooth NFW 
dark matter profile with a final annihilation state of bb 
in Fig. [3j we find that Coma yields poor limits for other 
dark matter profiles and annihilation channels as well. 

Finally, we searched for diffuse radio emission usin g the 
NRAO VLA Sky Survey (NVSS; ICondon et al1fl998h at 
1.4 GHz for some of the best candidate clusters for dark 
matter detection which do not have published radio up- 
per limits or detections, including Virgo, Fornax, AWM7, 
and two Virgo groups, M49 and NGC4636. We exclude 
Virgo and Fornax as they both contain very bright radio 
sources that likely wash out any true diffuse emission. 

For AWM7, we find no detected diffuse radio emis- 
sion in the NVSS images. We consider a circular emis- 
sion region with a diameter of ~ 2 core radii, and af- 
ter subtracting out point sources, we place an upper 
limit on the diffuse emission, following the method in 
Rudnick fc Lemmermanl (|2009l ). 
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Figure 2. Greyscale image of NGC4636 from NVSS, 45" resolu- 
tion, on a scale of 0-50 mJy/45" beam. Contours show the residual 
diffuse emission at a resolution of 135", as shown in the lower left, 
after subtraction of the compact emission. Contour levels are at 2 
mJy/135" beam x (1,2,3,...). 

We do find extended radio emission in the NVSS from 
M49 and NGC4636 extending beyond the bright emis- 
sion from the central galaxies, after removing the emis- 
sion from compact and slightly resolved sources using 
a mult iresolut ion filtering te chnique as used by e.g., 
iRudnick fc Lemmermanl (j2009l ). The resulting extended 
emission from M49 is approximately ~ 30 — 35 kpc 
(6 — 7') in diameter with a flux density of 40 mJy, much 
larger and fainter than the ~ 220 m Jy emission from 
the central, Seyfert 2 elliptical galaxy (jDunn et al.1 2010: 
iBrown et al.l l20llal) . NGC4636 shows diffuse emission 
over a ~ 36 kpc (7') range with a flux density of 30 mJy. 
The central source of NGC4636 show clear jet-like struc- 
ture and is approxim ately 5 kpc in size at 610 MHz 
iGiacintucci et all (j2011l ) , and this structure is not visible 
at 1.4 GHz at the lower resolution of the NVSS, whereas 
the diffuse emission at 1.4 GHz is larger and more spher- 
ical in shape. This diffuse central emission has not been 
previously identified in either group. Images of the dif- 
fuse emission for M49 and NGC4636 are shown in Figs. [I] 
and [21 respectively. 

3.2. Magnetic Fields in Clusters 

The observation of diffuse radio emission in galaxy 
clusters requires the presence of large-scale magnetic 
fields associated with the ICM. Intracluster magnetic 
fields can be inferred from the Faraday rotation of po- 
larized radiation of individual sour ces with typical cen- 
tral values of ~ 1 - 10 nG (e.g. lEilek fc Owenl l200l 
iBonafede et~aTI 120101 : IVacca et al.l I2012D . Clusters with 
cool cores can have higher inferred central ma gnetic 
field strengths, ~ 10 - 40 fiG (jTavlor et al.l I2002L 120061 : 
iKuchar fc Enf31inll2011f ). Magnetic fields in clusters can 
also be estimated using radio halo observations assuming 
equipartition of cosmic ray and magnetic energy densi- 
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ties; these estimates are typica lly between 0-1 — 1 /iG and 
are taken as lower limits (e .g. iGovoni Ferettill200l see 
also Carilli & Taylor 2002 for a review of magnetic fields 
in clusters). 

Turbulence due to past mergers can a mplify mag- 
netic fields in clusters to /iG levels (e.g.. iDolag et al.l 
2002: ISubramanian et al.ll2006URvu et al. 2008), and can 
also accelerate particles, which may be responsible for 
the giant radio haloes obs erved in some clusters (e.g., 
iBrunetti fc Lazari an 2011 and references therein). If 
mergers drive radio haloes, then the difference between 
clusters with and without haloes may be dynamical: 
clusters with haloes have suffered a merger in their re- 
cent past, while clusters without haloes are more dy- 
nami cally relaxed, with no recent mergers in their his- 
tory (IBrunetti et al,ll2QQ9b ICassano et al.ll2010D . Assum- 
ing clusters evolve in this way, it is natural to infer 
that at some point the cluster magnetic field will dis- 
sipate and suppress the sy nchrot ron radiation p roduc- 
ing haloes ([Brunetti et al.l 120091 : ICassano et "all 1201(1 ). 
However, theoretical studies and simulations of clus- 
ters show tha t large scale magnetic fiel ds are long-lived, 
~ 4 — 5 Gyr (Su bramanian et al.l 120061 ). which is longer 
than the cosmic ray electron/positron lifetime, and thus 
l onger than the lifeti me of the radio halo, about 1 Gyr 
(Bru netti et a l. 2009). Additi onal theoretical work on 
the ICM bv iKunz et al.l (j2011f h in which heating by lo- 
cal plasma instabilities can stabilize cooling in the ICM, 
predicts ~ 1 — 10 fiG magnetic field strengths that scale 
weakly with gas density and generally agree with values 
estimated from RMs for specific clusters. 

Only a handful of clusters have well-studied magnetic 
fields; less than half of our sample have published stud- 
ies on their magnetic fields. For the clusters that are 
not strong cool cores as defined by lHudson et al.l ((2010) 
but do not have any published information about their 
magnetic fi elds, we use the best fi t values of Bo and r] for 
Coma from Bona fede et al.l ([20101 ) , which we call the non- 
cool-core model, as Coma is a well-studied non-cool-core 
cluster. The magnetic field in clusters may turn out to be 
connected to the presence of a radio halo in such a way 
that clusters without observed haloes h ave lower mag- 
netic fields than those with haloes (e.g., IBrunetti et al.l 
2009). However, we choose to use the magnetic field pa- 
rameters derived for Coma even for clusters with upper 
limits on the radio emission. This is because RM studies 
yield no differences in the RMs, and thus inferred mag- 
netic fields, of clu s ters with and with out radio haloes 
(IClarke et al.ll200ll : IGovoni et al.ll2010l) . 

Cool core clusters generally have higher inferred central 
magnetic field strengths than clusters without cool cores 
(jKuchar fc Enfilinl 1201 ID . In our sample both A2029 
and Perseus have publishe d c entral field streng ths from 
RMs (jEilek fc Owenll2002l and lTavlor et al.ll2006l respec- 
tively). A2199 is the only cool core cluster in our sample 
with best fit values for both Bq and r] from Faraday RMs 
([Vacca et al.ll2012l ). Ophiuchus, also a cool core clus- 
ter, has a published RM from a radio galaxy outside the 
cool core region and therefore no estim ated central mag- 
netic field strength ([Govoni et al.l 120101). It has one of 
the hottest known cool cores (Fujita et al. 2008), which 
could mean its central magnetic field is somewhere be- 
tween that of a strong cool core cluster and a cluster with 
no cool core. We therefore choose the more conserva- 
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Figure 3. Dark matter annihilation cross section upper limits for 
all clusters in the sample; we assume a smooth NFW dark matter 
density profile. 

tive, non-cool-core model to estimate the magnetic field 
of Ophiuchus when comparing constraints from other 
clusters, but we also show how our constraints change 
if we use the central magnetic field value of Perseus, the 
prototypical cool-core cluster. We discuss the effects of 
the uncertainty on the magnetic field on dark matter 
constraints, specifically with respect to Ophiuchus and 
A2199, further in Sec. EH 

Little is known about the magnetic field strengths for 
groups of galaxies, especially in/near the central elliptical 
galaxy that typically dominates these groups. We choose 
to use the non-cool-core model for the magnetic field in 
calculating the limits for the two Virgo groups in our 
sample, M49 and NGC4636. These groups are heavily 
dominated by their central elliptical galaxies, which must 
host magnetic fields of their own. Large scale galactic 
magnetic fields of normal gala xies are typica lly in the fiG 
range (e.g., for the Milky Way. iNout sos 2012). We choose 
therefore to use the non-cool-core model for the magnetic 
fields of these groups, the central strength of which is 
probably comparable to or lower than the strength of 
the field within the central galaxy. Cluster properties, 
radio data, and magnetic field parameters are listed in 
Table El 

4. RESULTS AND DISCUSSION 

We derive constraints on the dark matter annihilation 
cross section by assuming conservatively that the radio 
upper limits or low levels of observed emission in a sample 
of nearby galaxy clusters are due to synchrotron emission 
from electrons and positrons produced by annihilating 
dark matter. In Fig. [3j we present the cross section up- 
per limits for all the clusters in our sample, each with an 
NFW dark matter profile, for the bb annihilation channel. 
For this dark matter profile, the Virgo groups NGC4636 
and M49 yield the tighest constraints, while the best 
constraints are produced by A2199 for massive clusters 
with published radio data at low dark matter particle 
masses and by Ophiuchus at higher masses. Perseus, 
which hosts a particularly bright mini-halo, produces the 
weakest constraints. As expected, the constraints from 
the Coma Cluster, with its bright, giant radio halo, are 
also relatively poor. 

_We also consider four different annihilation channels: 
66, r + r~, and In Figures [4] and [5] we 
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Figure 4. Dark matter annihilation cross section upper limits for 
A2199, for four annihilation channels, with a smooth NFW dark 
matter profile only. 

show the limits for each of the different channels for 
A2199 and Ophiuchus, using an NFW dark matter pro- 
file. The r + r~ and /i + \i~ channels, which tend to pro- 
duce more electrons and positrons per annihilation than 
bb, yield the best limits at lower particle masses, while 
the 66 yields better limits at higher masses, which is gen- 
erally true of all the clusters. 

4.1. Comparison to Limits from Gamma-Ray Emission 

Many previous studies of clusters have focused 
on the potential gamma-ray emission from dark 
matter annihilation, typically using data from 
Fermi jA ckerma nn et all l2Q10at lAndo fc Nagail [2Qll 
IHuang et all l2Q12t iNezri et al.l [20121 : lHan et all 120121 ). 
Using gamma rays to investigate the dark matter 
content of clusters has the advantage that the emission 
only depends on the underlying particle population(s) 
and chosen dark matter profile, and does not depend 
on other cluster properties, e.g., the magnetic field. 
However, only very nearby clusters produce useful 
limits, and clusters that are near the Galactic plane or 
contain a point source that is bright in gamma rays must 
be excluded from these studies, as the Galactic plane 
is gamma-ray bright and the resolution of Fermi and 
ground-based gamma-ray telescopes is poor compared to 
radio telescopes. A2199 and Ophiuchus, which yield the 
best limits in many cases in the radio, are typically not 
considered for gamma-ray analysis, as A2199 is too far 
away, and Ophiuchus is very near the Galactic center, a 
very bright gamma-ray region. Comparing the limits for 
A2199 for dark matter annihilating to 66 with a smooth 
NFW profile to the limit s derived for Fornax , also for 66 
and an NFW profile, in IHuang et al.l (|2012[ ). which are 
the best limits derived from gamma rays from a single 
cluster with an NFW profile to date, our limits from 
A2199 are approximately a factor of ~ 3 or more across 
a wide range of m asses. 

Han et al. (2012), using; the results of the Phoenix sim- 
ulations to model substructure and considering the pos- 
sibility that the gamma-ray emission from clusters is ex- 
tended, placed constraints on dark matter annihilation in 
Coma, Virgo and Fornax t hat are compa r able t o or better 
than the limits derived in IHuang et all (|2012[ ) including 
substructure. Our limits from A2199 when considering 
the same substructure model from the Phoenix simu- 
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Figure 6. Dark matter annihilation cross section upper limits 
for A2199, for one annihilation channel (bb) and our four different 
substructure models. 



lation s, are comparable to those reported in lHan et all 

4.2. Substructure in Clusters and Groups 

The boost factors obtained by using the the Phoenix 
simulations are > 10 times larger than the boost fac- 
tors derived from our other adopted substructure models, 
which in turn yield limits that are 1-2 orders of magni- 
tude tighter. This is reflected in Fig. [6j which shows the 
limits for the bb annihilation channel for our four adopted 
dark matter profiles for A2199. The limits from our con- 
servative and optimistic models are only slightly better 
than those produced by an NFW profile, about 10 — 30% 
. The limits that result from using the Phoenix profile 
are much lower, almost 2 orders of magnitude for A2199, 
and dip below the nominal thermal annihilation cross 
section 3 x 10 -26 cm 3 s _1 for masses < 400 GeV. A sim- 
ilar pattern is true for most other clusters, except M49 
and NGC4636. This is because the emission region is so 
small for these groups that the J factors only change by 
at most a factor of four as the amount of substructure is 
increased, while for other clusters the difference between 
J nfw and Jphx is typically 2 — 3 orders of magntitude. 
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Figure 7. Effects of uncertainty in cluster mass and magnetic 
field parameters, for A2199, for the final state 66, with an NFW 
dark matter profile. The constraints from the Fornax Cluster using 
gamma ray u pper limits from Fermi are copied from Fig. 4 of 
Huan g~et al.l (| 20 121 ) and also correspond to an NFW profile for 
final state 66. 
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Figure 8. Effects of uncertainty in the central magnetic field 
strength for Ophiuchus, for one annihilation channel (66), with an 
NFW dark matter profile. The constraints from the Fornax Cluster 
us ing gamma ray upp er limits from Fermi are copied from Fig. 4 
of Huang et al. ( 2012) and also correspond to an NFW profile for 
final state 66 (same as in Fig. 0. 



4.3. Uncertainties in Cluster Masses and Magnetic 

Fields 

Aside from the assumed dark matter profile and 
amount of substructure, the two main sources of uncer- 
tainty in the limits on the dark matter annih ilation cross 
section are the cluster mass, reported by iChen et al.l 
(|2QQ7[ ) , and the uncertainty in the magnetic field param- 
eters Bq and r\. IChen et all (j2007f ) also report uncertain- 
ties in the gas fraction, the core radius, and /?; we find 
that the uncertainties associated to these quantities are, 
however, negligible compared to the cluster mass and 
magnetic field uncertainties. 

In Fig. [7J we show the relative magnitudes of the var- 
ious sources of uncertainty for one of the best cluster 
candidates, A2199. The uncertainty in the mass leads to 
an uncertainty in the J factors which is generally < 2 
smaller for all the clusters in our sample. For the uncer- 
tainty in the magnetic field of A2199, we choose values 
for Bo and n that best rep resent the sp r ead in the Bq-t] 
contour plot in Figure 9 of iVacca et all {2012). The un- 
certainty in the magnetic field is comparable to the un- 
certainty in the cluster mass at lower dark matter masses, 
leading to uncertainties in the annihilation cross section 
of about a factor of 2. At higher masses, the magnetic 
field uncertainty translates to only about a 10% uncer- 
tainty in the limits while the uncertainty in cluster mass 
remains constant. This is because the total radio energy 
per frequency, and thus the annihilation cross section, 
depends on the magnetic field to a power proportional 
to the index of the underlying electron/positron distri- 
bution. The electron/positron distribution produced by 
lower mass dark matter particles is generally steeper than 
the distribution produced by higher mass dark matter 
particles, so the magnetic field is more important at lower 
dark matter masses than at higher masses. 

A2199 is a cool-core cluster (e.g.. iHudson et aDl2010f ) 
and is therefore more likely to host a small mini-halo hun- 
dreds of kpc in size rather than a giant, Mpc-sized radio 
halo, as we assumed h ere based on the upper limits a vail- 
able in the literature (Rudnic k fc Lemmermanll2009[ ). It 
is also poss i ble th at the magnetic field determined by 
IVacca et al.l (j2Q12h is only relevant inside the cool core 



region. However, a smaller emission region would likely 
not change our results much, since the dark matter distri- 
bution and thus emission from dark matter is highly cen- 
trally peaked (except perhaps in the case of the Phoenix 
simulations, where the dark matter distribution is much 
flatter). 

While Ophiuchus yields some of the best limits for 
the r + r~ and ji + \i~ annihilation channels, its mag- 
netic field is not well understood. Ophiuchus has a 
claimed detection of nonthermal ha rd X-ray emission 
from a d eep INTEGRAL observation (jEckert et al.ll2008l : 
see also iProfumol l2008f ) . The magnetic field strength 
of Ophiuchus has been estimated to be around 0.1 fiG 
using this observation combined with radio data, as- 
suming the hard X-ray emission is due to the IC scat- 
tering of the same population of relativistic electrons 
that are the sou rce of the radio emission due to syn- 
chrotron losses ( Eckert e t al.l 120081 : iPerez-Torres et all 
2009: iNevalainen et a l. 2009). This is inconsistent with 
the larger magnetic fields inferred from Faraday RMs in 
similar clusters. 

Both approaches to measuring cluster magnetic fields 
rely on several assumptions. The field strength calcu- 
lated using hard X-ray emission is typically a volume av- 
erage over the size of the emission region, and it assumes 
that the field is uniform, while simulations point to a field 
that decays radially with gas density (e.g. JPonnert et all 
2009). Additionally, assumptions must be made about 
the energy spectrum and origin of the cosmic ray lepton 
population in order to match the observed X-ray and ra- 
dio emission that potential ly conflict with other observa- 
tions , e.g. in gamma rays (Colafr ancesco fc Marchegianil 
[20091). 

RMs rely on the gas density along the line of sight. 
Turbulence could distort the gas density and therefore 
the magnetic field along the line of sight. Additionally, it 
is possible that most of the contribution to an observed 
RM is due to the local environment around the radio 
source itself, implying the inferred magnetic fields may 
not be representative of w hat is happenin g in the ICM 
dRudnick fc B lundell 2003, however see also lEnfilin et al.l 
2003]). 
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We choose to use field estimates based on RMs of other 
clusters for Ophiuchus, setting the non-cool-core model 
as the lower limit, and the field of Perseus, a cluster with 
a strong cool core, as the upper limit. We show the effect 
of the uncertainty in the magnetic field of Ophiuchus in 
Fig.[8]by varying Bo with r\ fixed at 0.5. The limits vary 
by approximately a factor of 7 for a mass of 10 GeV, and 
a factor of 1.3 at 1000 GeV. 

5. CONCLUSIONS 

We presented here new limits on the dark matter an- 
nihilation in clusters of galaxies using radio data in clus- 
ters. Overall, the limits derived from upper limits on 
radio haloes in massive clusters, low levels of observed 
central diffuse emission in galaxy groups, and detections 
of mini-haloes in cool core clusters are better than previ- 
ously derived limits using non-detections in gamma rays 
for the same substructure model and annihilation chan- 
nel by a factor of ~ 3 or more. However, our limits de- 
pend strongly on both the assumed amount of substruc- 
ture, which is also true for limits derived from gamma 
rays, and the chosen magnetic field model for any in- 
dividual cluster. We consider four models for the dark 
matter spatial profile, with an increasing level of sub- 
structure. Using the results from the Phoenix simula- 
tions, which contain the highest fraction of substructure, 
the limits on the annihilation cross section are 1—2 orders 
of magnitude better than those adopting an NFW profile, 
and lower than the thermal relic cross section for masses 
< 400 GeV. The uncertainties in the magnetic field mod- 
els for A2199 and Ophiuchus, which yield the best lim- 
its for a range of annihilation channels and substructure 
models, lead to uncertainties in the cross section limits 
of a factor of 2 — 7 for low mass dark matter particles and 
a factor of 2 or less for higher masses. We consider four 
dark matter pair-annihilation channels in our analysis; 
channels that produce harder leptons (r + r~ and 
yield better limits at lower masses, as expected, while at 
higher masses, bb yields the best limits. We also report 
the detection of faint diffuse radio emission from M49 
and NGC4636, two groups near the Virgo Cluster, using 
NVSS data. These groups yield very strong limits for an 
NFW dark matter profile, but since the emission regions 
are so small, they only improve by a factor of < 2 as the 
amount of substructure is increased. 
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is supported in part by NSF grants AST-0908688 and 
AST-1211595. 
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Table 1 
Galaxy Cluster Properties 



Name 




z 






P 


r c 


V 


S u 




Radio 


Bo 


Bo 








(10 14 M ) 






(kpc) 


(MHz) 


(Jy) 


(Mpc) 


Ref. 


(/iG) 


Ref. 


A0576 





0381 


4.80 





825 


277 


327 


<0.20 


0.5 


1 






A2199 





0302 


4.18 





655 


98 


327 


<0.025 


0.5 


1 


11.7 


4 


A2244 





0970 


5.14 





607 


88 


327 


<0.08 


0.5 


1 






Zwl742 a 





0757 


11.1 





717 


163 


327 


<0.10 


0.5 


1 






Perseus 





0183 


5.04 





540 


44 


1400 


1.979 


0.069 


2 


25 


5 


Ophiuchus 





0280 


48.4 





747 


196 


1400 


0.1064 


0.315 


2 






A2029 





0767 


9.62 





582 


58 


1400 


0.0188 


0.125 


2 


16.0 


6 


Coma 





0232 


9.88 





654 


241 


1400 


0.64 


0.415 


3 


4.7 


7 


AWM7 





0172 


4.92 





671 


122 


1400 


<0.107 


0.122 


* 






M49 




b 


0.67 





592 


7 


1400 


0.040 


0.016 


* 






NGC4636 




c 


0.18 





491 


4 


1400 


0.030 


0.018 


* 







Notes. 

a Full name for Zwl742 is ZwC11742. 1+3306. 

b Distance calculated from a distance modulus of 31.15, from Villegas et al. (2010), rather than a redshift. 
c Distance calculated from a distance m odulu s of 31.24, from Dirsch et al. (2005), rather than a redshift. 

Cluster masses taken fro m Chen et al. (2007), corrected for gas fraction and current cosmology (see text). Listed values for z, (5 and r c 
also taken from Chen et al. (2007) (r c corrected for current cosmology). 

For all clusters except A2199, we assume rj = 0.5. We use 77 = 0.9 for A2199, consistent with IVacca etaTI (2012). 
For clusters with no listed Bo, we assume Bo = 4.7/xG and 77 = 0.5, i.e., our non-cool-core model. 
References. 

* Thi s pa per. (1) |R udnick & Lemmerman (2009) (2) Mur gia et all (|2009h (3) IKim et al.l jl990h (4) IVacca et all (|20ll ) (5) ITavlor et all 
( 2006) (6) Eilek & Owen ( 20021) (7) IBonafede et all (|2010l ) 



